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(Dated:) 

We analyze the critical connectivity of systems of penetrable d-dimensional spheres having size 
distributions in terms of weighed random geometrical graphs, in which vertex coordinates correspond 
to random positions of the sphere centers and edges are formed between any two overlapping spheres. 

Edge weights naturally arise from the different radii of two overlapping spheres. For the case in which 
the spheres have bounded size distributions, we show that clusters of connected spheres are tree-like 
for d —> oo and they contain no closed loops. In this case we find that the mean cluster size diverges 
at the percolation threshold density r/ c —> 2 ~ d , independently of the particular size distribution. We 
also show that the mean number of overlaps for a particle at criticality z c is smaller than unity, 
while z c —y 1 only for spheres with fixed radii. We explain these features by showing that in the 
large dimensionality limit the critical connectivity is dominated by the spheres with the largest size. 

Assuming that closed loops can be neglected also for unbounded radii distributions, we find that the 
asymptotic critical threshold for systems of spheres with radii following a lognormal distribution is 
no longer universal, and that it can be smaller than 2 -d for d —» oo. 


I. INTRODUCTION 


Percolation phenomena are ubiquitous in many aspects 
of natural, technological, and social sciences, and they 
arise when system-spanning clusters or components of, 
in some sense, connected objects form [II- A quan¬ 
tity of much interest is the percolation threshold, which 
marks the transition between the phase in which a gi¬ 
ant component exists and the one in which it does not. 
In general, the percolation threshold is a nonuniversal 
quantity, as it depends on the connectivity properties of 
the specific system under consideration Q. For example, 
in continuum percolation systems, where objects occupy 
positions in a continuous space, the threshold depends 
on the shape of the objects [Jt 3, on their interactions 
M , as well as on the connectedness criteria EMU- 

In this article, we consider the infinite-dimensional 
limit of a paradigmatic example of continuum perco¬ 
lation: the Boolean-Poisson model [3, HU- In this 
model, penetrable spheres with distributed radii have 
centers generated by a point Poisson process, and any 
two spheres are considered connected if they overlap. For 
a given distribution of the radii, the percolation thresh¬ 
old is given by the critical concentration r) c of spheres, or 
by the critical volume fraction (f) c = 1 — e ~ Vc , such that 
a giant component of connected spheres first forms. Pre¬ 
cise numerical estimates of r) c have been obtained in two 
and three dimensions for, respectively, disks and spheres 
with fixed or distributed radii [l6l - l23| . The general trend 
observed by these investigations is that r] c depends on the 
form of the distribution function of the radii, and that it 
has its minimum when the sphere radii are monodisperse 
(i.e., when the spheres have identical size)^ This last 
point has been formally confirmed in Ref. [3]) although 
it may not hold true in the limit of infinite dimensions d 

mm 

Here we show that for bounded distributions of the 
radii, that is for polydisperse spheres with a maximum 
finite value of the radius, the percolation threshold of the 


Boolean-Poisson model tends asymptotically to a univer¬ 
sal constant as d —> oo, provided that the radii distribu¬ 
tion is independent of d. This constant coincides with 
the value found in Refs. j.27i, !28j for spheres of identical 
radii, r] c —> 2~ d , and it is independent of the particu¬ 
lar form of the size distribution function. We interpret 
the universality of ?y c as being due to the statistical ir¬ 
relevance of the spheres with smaller radii: the onset of 
percolation is established effectively only by the subset 
of spheres with maximum radius. Furthermore, we show 
that the mean number of connected spheres per particle 
at percolation, z c , is less than unity for polydisperse dis¬ 
tributions of the radii, while z c —> 1 only in the limit of 
identical radii. This finding is analogous to what simu¬ 
lations have shown for the case of continuum percolation 
in three dimensional space of spherocylinders with length 
polydispersity |29j . 

These results rest on the observation that closed loops 
of connected spheres can be neglected in the limit of large 
dimensions, as we show explicitly for the case of bounded 
radii distributions. In the hypothesis that closed loops 
are irrelevant also for spheres of unbounded size, we show 
that 7? c for d —i oo is not universal, as it depends on the 
parameters of the distribution, and that it can be smaller 
than the critical threshold of monodisperse spheres, in 
contrast to what is expected for finite dimensions [24| . 


II. THE MODEL 


To construct the Boolean model, we consider N points 
placed independently and uniformly at random in a d- 
dimensional volume V. Each point is the center of a 
sphere with the radius drawn independently and ran¬ 
domly from a given probability distribution function 
p{R). If we denote N\ the number of spheres of radius 
J?i, N% the number of spheres of radius R 2 , and so on, 
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FIG. 1 . (Color online) Connectedness criterion for spheres 
with different radii, (a) The spheres of radii Ri and R 2 
overlap the sphere of radius R 3 , forming links between R 3 
and Tii and between R 3 and T? 2 . (b) Corresponding cluster 
formed by nodes (sphere centers) labeled by the sphere radii 
and weighted links (solid lines) connecting the nodes. 


we can write the following without loss of generality: 

p(R)=Y / xA R ~ R i), W 

i 

where Xi = Ni/N with i = 1,2,... is the fraction of 
spheres of radius Ri. 

Given any two spheres of radii, say, Ri and Rj , we 
assign a link between their centers if the spheres overlap, 
that is, if the distance r between their center is smaller 
than Ri+Rj, as shown in Fig. [I] We express this criterion 
for the formation of a link in terms of the connectedness 
function: 


f ij (r) = 6(R i +R j -r), (2) 

where 6 (x) = 1 for x > 0 and 9(x) = 0 for x < 0 is the 
unit step function. 

The set of sphere centers (nodes) and links (edges) 
forms a type of weighted random geometric graph, in 
which the probability that an edge between two nodes is 
formed is weighted by the sphere radii. To see this, let 
us take a sphere of radius Ri centered at the origin. The 
probability that a second sphere of radius Rj forms a link 
with the first sphere is: 

<i = yj = n d {Ri+ v Rj)d , ( 3 ) 

where dr is an infinitesimal d-dimensional volume el¬ 
ement at the position r of the sphere of radius Rj, 
f Id = 7T d / 2 /r(l + d/ 2 ) is the volume of a sphere of unit 
radius, and T is the gamma function. We note that 
defines also the excluded volume 1W = fld(Ri + Rj) d in 
units of V between two spheres of different radii. 


III. IRRELEVANCE OF CLOSED LOOPS FOR 

d —» 00 

An important aspect of the topology of random geo¬ 
metric graphs is represented by closed loops (or cycles) 


of connected nodes. The most studied loop quantity is 
the three-nodes cycle c d , often denoted the cluster coef¬ 
ficient, which gives the conditional probability that two 
nodes are connected given that both nodes are connected 
to a third one. c d has been calculated for systems of 
spheres with identical radii and for any dimensionality 
U The observation that c d ^ vanishes exponen¬ 

tially as d —> 00 indicates that random geometric graphs 
in large dimensions have a locally tree-like structure. 

Using results from the theory of hard-sphere fluids, it 
is actually possible to show that, in the limit of large di¬ 
mensions, closed loops are negligible also for any number 
of nodes and for bounded radii distributions. Random 
and weighted random geometric graphs have thus tree¬ 
like structures when d —> 00 . To see this, let us first 
consider the case of monodisperse spheres with radius 
Rm ■ We define an n-chain graph as a cluster of n > 3 
nodes with n — 1 edges such that every two consecutive 
edges, and only those, have a common node. We denote 
as end-nodes the two nodes of an n-chain that each have 
only one edge. The n-cycle coefficient c^ is defined as 
the conditional probability that two nodes are connected 
given that they are the end-nodes of an n-chain. Since 
the spheres have identical radii, we omit the subscripts 
in Eq. m, and we write the connectedness function as 
simply f(r) = 9{2 Rm — r). From the definition of c" d \ 
we can thus write: 

(n) = Jdr( n) f {|ri -r 2 |)/(|r 2 -r 3 |)---/(|r ra -ri|) 

/ dr(")/(|ri - r 2 |)/(|r 2 - r 3 |) • • • /(|r„_i - r n |) ’ 

(4) 

where dr^ = dridr 2 ■ • • dr n . Besides a prefactor, the 
above expression coincides with the cluster integral of 
a ring of n hard-spheres of radius Rm [SI], as the 
Mayer function JuM for a fluid of hard-spheres is just 
/mM = —fi r ) 0, IQ. To evaluate Eq. © for d —> 00 , 
we thus use known results from the theory of hard-sphere 
fluids in infinite dimensions. Noting that the denomina¬ 
tor of E q. (HI) (i.e., the n-chain contribution) is simply 
VV ^- 1 [33}, where V ex = VLd‘2- d R < M i s excluded vol¬ 
ume for spheres of identical radius Rm- and introducing 
the Fourier transform /(q) of the connectedness function 
we rewrite Eq. 0 as: 



yn-l 

Fex 


dq 

( 2 n) d 


/(q)”- 


(5) 


The integration in Eq. © for d —> 00 has been worked 
out in Ref. 33 ( see a l so Ref- H), so that the n-cycle 
coefficient reduces to: 


„(«) 


/ n — 2 

( n ) 

n/2 

n"' 2 

T—1 

1 

\n- 2 ) 

_(n - l) ra-1 _ 


-1 d /2 


, ( 6 ) 


from which we see that closed loops of any number n 
of nodes are exponentially small as d —> 00 , because the 
quantity within square brackets is less than unity for n > 
3. 
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Let us now consider the n-cycle coefficient for the case 
of polydisperse spheres. Using Eq. @ for the connected¬ 
ness function, we generalize Eq. 0 as follows: 




(™) _ 


(cfa 

Kit 


, l n 


,l n 



( 7 ) 


where 


c iu...,i n = J dr (n) /iij 2 (|ri -r 2 |)/ i2i3 (|r 2 -r 3 |)--- 




• X /inii(l r n - r l|), 

(8) 

v (n) . 

= / dr ( ” ) / ili2 (|ri-r 2 |)/ I2i3 (|r 2 -r 3 |)- 




• X /in-li„(l r »-l - r„|), 

(9) 

and 




((• 

))ii,...,i„ — 

F, X il X i 2 ---Xj n (---) 

(10) 



il in 


denotes 

a multiple 

average over the 

radii 


.,R in . In 

the appendix, we show 

that 


for bounded distributions of radii, the n-cycle coefficient 
in the limit d —> oo is such that: 

(Q) W < c^xP, (11) 

where is the n-cycle coefficient for identical radii, 
Eq. db]), and Xd ' 1 oc d a , where a is a nonnegative constant. 
Since the exponential decay of c ^ n ' 1 for d —> oo is stronger 
than the power-law increase of ^ , we see thus that 
also for the case of polydisperse spheres for bounded radii 
distributions, the n-cycle coefficient vanishes for any n > 
3. 


IV. SIZE OF FINITE COMPONENTS 

The observation made in the previous section that 
closed loops are irrelevant in the large dimensional limit 
of the Boolean model allows us to consider the compo¬ 
nents of the associated weighted random geometric graph 
as effectively having a tree-like structure. This leads to a 
considerable simplification, as we can take the formalism 
of the theory of random graphs (see, e.g., Refs. [35l - l37| ) 
and generalize it to the case in which nodes have weights. 


A. Multidegree distributions 

We start by considering the multidegree distribu¬ 
tion of a node of type i, defined as the probability 
Pi(l, fa-, 2, k 2 ;...) that a sphere of radius Ri is connected 
to fci spheres of radius R\, k 2 spheres of radius i? 2 , 
and so on. Since the radii are randomly and indepen¬ 
dently distributed among the N nodes, Pi(l, fa-, 2, fc 2 ;...) 


is just a product of binomial distributions pij(kj) (with 
j = 1,2,...), each giving the probability that kj spheres 
of radius Rj overlap the sphere of radius Rp. 

Pi{l, fa-,2, k 2 -,...) = Y[pij{fa), (12) 

j 

with 

P*j(k J )=( Nj -' Si ’ j yvli) k fa l-v£) N *- s '-*~ k *, (13) 

where Nj (with j = 1,2,...) is the number of spheres 
of radius Rj, vl J x are the overlap probabilities given in 
Eq. ([3]), and Sij is the Kronecker symbol. 

We next consider for all i the limit Nj —> oo such that 
Nj/V — Xip remains finite, where p = N/V is the to¬ 
tal number density. In this limit, Eq. Cl reduces to a 
Poisson distribution: 

kj 

Pij{kj) = ~j7J e ( 14 ) 

where 

Zjj = ^2 kpij(k) = Xjpn d (Ri + Rj) d (15) 

k 

is the average number of spheres with radius Rj that 
overlap a given sphere of radius Rj. 

In addition to the node degree distribution 
Pi(l, fa-,2, k 2 ', ■■ ■), for the following analysis we 
will also need the excess node degree distribution 
qjiil, fa-,2, k 2 -,...), defined as the conditional probability 
that a sphere of radius Rj is connected to fa spheres of 
radius Ri (with l = 1,2,...), given that it is connected 
to a sphere of radius Rj. This task is simplified by the 
irrelevance of closed loops in the large dimensionality 
limit. In this case, indeed, if we select at random an 
edge connecting a node of type j with a node of type 
i, the j node attached to the edge is fa times more 
likely to have degree fa than degree 1 with nodes of type 
i. Its degree distribution will thus be proportional to 
fa/Pj ( 1 , fa-,2, k 2 ; ■■ ■)■ The excess degree distribution of a 
j node that has fa edges with nodes of type i other than 
the edge with the node i to which is attached is thus 
0: 


qji{l, fa-, 2 ,k 2 , ■ 

where = 

qji{l, fa-, 2 ,k 2 ,. 


_ (fa + l)pj(l, fai, fa + 1 ;...) 

Ek( k i+ 1 )Pj( 1 . fc i; •••;*. fe+i; ■ • ’ 

(16) 

Efci, fc2 ,...- From Eqs. C3 and (Jl4j), 

..) reduces simply to: 


Qji{ 1? 2, • • •) 


(fr+ 1 )^ + 1 ) ^ pAkl) 

Zji , ,. 

Y[pji(kt), ( 17 ) 


where we have used (fa + l)pji(fa + l) = ZjiPji(fa). Equa¬ 
tion (ED states thus the well-known property that the 
excess degree distribution coincides with the node degree 
distribution when this is Poissonian [35j . 
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FIG. 2. (Color online) Schematic representation of a finite 
tree-like cluster formed by connected spheres of radii I?i, R 2 , 
and R 3 . Each node label corresponds to the value of the 
radius of the sphere attached to the node. 


B. Mean cluster size in the subcritical regime 

We exploit now the statistical irrelevance of closed 
loops discussed in Sec. 1331 ] to find the mean size S of 
finite clusters of connected spheres as d —> 00 . In doing 
so, we shall first keep the form of the degree distributions 
unspecified, and apply Eqs. © and © only at the end 
of the calculation. 

Let us start by considering a randomly selected node 
that has probability Xi of being occupied by a sphere of 
radius Ri. Due to the general tree-like structure of the 
graph, the cluster to which the selected node belongs is 
formed by branches attached to the node according to 
the degree distribution pt(l, k\\ 2, fc 2 ; • • ■), as schemati¬ 
cally shown in Fig. [2] The mean size Si of the cluster to 
which the selected node belongs is thus: 

Si — Xi T Xi ^ ( Pi (1, k\ , 2, A: 2 , • * •) ^ [ kj Tij , (18) 

k j 

where T.\j is the mean cluster size of one of the kj 
branches attached to the selected node. Since the clus¬ 
ters have a tree-like structure, T.\j is given by the mass 
(unity) of one neighbor of the selected node, plus the 
mean cluster size of each of the remaining subbranches 
attached to the neighbor. To find Rj. we thus need the 
excess degree distribution qji(l, fci; 2, fc 2 ; • • ■) of a sphere 
of radius Rj connected to the selected node of type i: 

= 1 + ^ qjj( 1, fci; 2, fc 2 ; ■ ■ ■) y"] kjTji. (19) 
k 1 

Equations (fl8l) and fijil) are quite general, as they ap¬ 
ply also to tree-like graphs with degree distributions that 


are not reducible to a multiplication of Poissonian prob- 
abilitites. Interestingly, similar equations are found in 
the calculation of finite size components of multigraphs 
(also denoted multiplex networks), formed by different 
networks, each having particular node properties, cou¬ 
pled together |H, [3§J. The Boolean-Poisson model with 
random radii can thus be viewed also as a particular type 
of multigraph, in which each individual network is con¬ 
stituted by nodes occupied by spheres of a given radius. 

Let us now use the results of Sec. IIV Al and rewrite 
Eqs. m and m by substituting pi( 1, k\\ 2, /c 2 ;...) and 
qji( 1, fci; 2, fc 2 ;...) with, respectively, Eqs. (fl2l) and (flTl) : 

Si = Xi + Xi ^2 E (k)Tj = Xi + Xi ^2 z ijTj, (20) 

3 k 3 

= 1 + E E k P3i(k)T = 1 + J2 ZjiTu (21) 

l k l 

where we have used Eq. m and the fact that Rj de¬ 
pends only on the neighbor ( j) of the selected node, i.e., 
T =T 

The mean cluster size is given by S = which 

from Eq. (l20l) reduces to: S = 1+£T . XiZijTj. This rela¬ 
tion is obtained also if we multiply both sides of Eq. m 
by Xj and sum over j. We can thus write: 

S = 22 x i T i' ( 22 ) 

3 

which states that S is just the average over the sphere 
radii of the mean cluster size of the branches. 


C. Equivalence with the Ornstein-Zernike equation 
for the pair-connectedness 

In continuum percolation theory, cluster statistics are 
often studied using the formalism of pair-connectedness 
correlation functions [1, [51 HjJ , which exploits well devel¬ 
oped techniques of liquid state theory. This method has 
been recently used to studying percolation of monodis- 
perse spheres in large dimensions [13 • 

As long as closed loops can be neglected, the network 
formalism discussed above and the pair-connectedness 
functions method give identical results, provided that the 
second-virial approximation is taken. To see how this 
equivalence holds true for the Boolean model in large 
dimensions, let us first consider the pair-connectedness 
function Pij{r — r'), defined such that XiXjp 2 Pij( r — 
v')dvdv' is the probability of finding two spheres of radii 
Ri and Rj within the volume elements dr and dr' cen¬ 
tered respectively in r and r', given that they belong to 
the same cluster. The mean cluster size S is given in 
terms of Pij(r — r') by the following relation[42j: 

5 = 1 + p^2 XiXjPij, 

i,3 


( 23 ) 
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where Pjj = f drPjj (r). Pjj is the solution of the pair 
connectedness analog of the Ornstein-Zernike equation of 
the liquid state theory of fluids: 

P%j = Cij + p y ' XiCgPij, (24) 

i 

where Cij = J drCij(r) is the volume integral of the di¬ 
rect pair connectedness function C'y(r), which describes 
short-range connectivity correlations. Let us introduce 
the quantity T) defined as: 


valid. To measure the sphere concentration we introduce 
the dimensionless density 

V = P n d (R d ) R = P n d XiRf. (29) 

i 

The percolation threshold rj c is defined as the smallest 
value of 77 such that S diverges. This definition is equiv¬ 
alent to finding the smallest pole of Eq. ED, if it exists. 

A. Discrete radii distributions 


Ti = 1 + p'Y^x j Pij. 


(25) 


The use of the above expression reduces Eq. (1231) to: 


S — ^ ( Xi T p ^ ) XiXjPij — ^ ( Xi I 1 T p ^ ( XjPij 
i i,j i \ j 

= ^2,Xifi, (26) 

while inserting Eg. (1211) into Eq. (1251) leads to: 

P = 1 + p^Xj ^'Cij + p^XiCuPij^j 

= l + p£ XjCij + P£ XjCij C Tj - 1) 

3 3 

= 1 + pY / x j C ij T j . (27) 


We first consider the case in which the spheres have a 
finite number M of radii: 

M 

p{R) = Y,XiS(R-Ri), (30) 

so that using Eqs. m , ED and (1221) we rewrite the 
equations for the mean cluster size as: 

M 

S = J2xi T i, (31) 

2=1 

M 

T t = l+pn d J2 Xj {Ri + Rj)%. (32) 

j =1 

Without loss of generality, we assume that Rm is strictly 
the largest radius out of the M possible values of the 
radii, and we introduce qi = Ri/R m , which takes values 
smaller than the unity for all i yf M. For large d , the 
dimensionless density 77 reduces to: 


We see that Eqs. ESD and ED are identical to respec¬ 
tively Eqs. (1221) and ED if we identify pXjCij with Zij. 
From Eq. CD, we obtain thus: 

Cij = = n d (Ri + Rj) d , (28) 

pXj 

which corresponds to take the volume integral of 
the second-virial approximation Cy(r) = Gf'(r) = 
fij{r) for the direct pair-connectedness function. This 
is not surprising, because in the density expansion 
of the direct pair-connectedness function, CV,(r) = 

J2 n >2 p" _ 2 C/(£r), the terms with n > 3 contain at least 
one closed loop. 


V. UNIVERSALITY OF THE PERCOLATION 
THRESHOLD 


M 

77 = pfl d y XiR d 
■A pn d RMXM, 


P^clRm 


M—l 

XM + £ XiQ d 
i=1 


(33) 


because q d goes exponentially to zero as d —> 00 when 
i 7 Z M, and Eq. (1521) becomes: 


Ti = 1 + 2 d 77— y xj 
XM “ 

3 


{H 31 ) T ‘ 


(34) 


—> 00 


We note that ^ i s vanishingly small as d 

unless i = j = M , for which it takes the value 1. The 
smallest pole of Eq. (1511) for large d is thus the solution 
of: 


Ti = 1 + 2 d v T M 6i, M , (35) 


We proceed to find the percolation threshold for the 
Boolean-Poisson model of polydisperse spheres in the 
large dimensionality limit. We shall consider the case 
of bounded distributions of the radii, for which we have 
shown in Sec. mu that closed loops of connected particles 
can be neglected for d —> 00 , and Eqs. ED and (1221) are 


where Sjj is the Kronecker delta. Equation (1351) is solved 
by T m = 1/(1 — 2 d rj) and T) = 1 for i ^ M, so that the 
mean cluster size ED becomes: 

M—l 

S=y Xi + x m Tm = l ' VM 2 d rj , (36) 
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FIG. 3. (Color online) (a) Percolation threshold r/ c in units of 
the asymptotic value 2 ~ d as a function of dimensionality for 
a discrete distribution of radii with M = 2 and R 1 /R 2 = 1/2. 
In this case Eq. (152)1 is a system of two linear equations so 
that r/c can be calculated exactly for any d. X 2 = 0.8, 0.6, 0.4, 
and 0.2 from the uppermost to the lowermost curves, (b) 
Corresponding values of the critical coordination number z c . 
As d —> 00 , z c tends asymptotically to x 2 - 


which diverges when 


1 

V^rj c =^. (37) 

The above expression for rp. holds true for any sequence 
of occupation fractions Xi , independent of dimensional¬ 
ity, provided that % ^ 0. In particular, Eq. 43ZD con¬ 
firms and extends to M > 2 the finding of a previous 
report that spheres with two different radii have a uni¬ 
versal critical threshold in infinite dimensions [26| . Note 
that r/c = 1 /2 d is also the limit for infinite dimensions of 
the percolation threshold of monodisperse spheres with 
radius Rm , whose mean cluster size is given by Eq. 
with xm = I- 

The origin of the universality of r] c can be traced back 
to the divergence of Tm, which indicates that the onset of 
a giant component of connected polydisperse spheres is 
established only by the subset of spheres with the maxi¬ 
mum radius when cL —> 00 . In other words, at d —> 00 the 
contribution to percolation of the smaller spheres van¬ 
ishes, and the resulting r] c is the critical threshold for a 
system of monodisperse spheres of radius Rm- Follow¬ 
ing the observation that systems of polydisperse spheres 
with M different radii can be interpreted as a multinet¬ 
work of coupled M subnetworks (see Sec. HVBl) , we see 
that Eq. (|35l) is equivalent to decoupling the subnetworks 
associated with each radius, and that long-range connec¬ 
tivity arises only from the network formed by spheres 
with radius Rm- 

One interesting consequence of the irrelevance of 
smaller radii at percolation is that the critical average 
connectivity per particle, 

z c = '^2 XiXjp c <d d (Ri + Rj) d , (38) 

h3 


reduces for d —> oo to: 

Zc — 2 ^ ^ x i x j f — 2 ~ j ^ ^ Vc x M = X M 5 

i, 7 ' ' 

(39) 

where p c is the critical number density. For xm < 1, the 
critical average connectivity is thus less than the unity 
for d —> 00 , which must be contrasted to z c > 1 for 
systems constituting only of monodisperse spheres in any 
dimension [28j . 

For the binary case (M = 2), Eq. (1521) reduces to a 
system of two linear equations that can be solved exactly 
for any d. The resulting ry c and z c are shown in Figs.[3Ka) 
and [3Kb), respectively, for R 2 = 2i?i and different values 
of the fraction x 2 of spheres of radius R 2 . The asymptotic 
limits T] c = 2~ d and z c = x 2 are recovered for sufficiently 
large values of d. 


B. Continuous radii distributions 


Let us now consider the case in which the radii distri¬ 
bution p(R) is a continuous bounded function indepen¬ 
dent of d. We again denote by Rm < 00 the maximum 
allowed radius, so that p(R) = 0 for R > Rm, and we 
rewrite the equations for the mean cluster size in terms 
of continuous variables of the radii: 


5 = {T(R)) r , (40) 

T(R) = l + P n d {(R + R') d T(R')) R :, (41) 

where ((• • • ))u = JdRp(R)(- ■ ■ ). We expand the bi¬ 
nomial power (R + R') d and use rj = pVLd{R d )n to write: 


T(R) = l+,^Q^A(R‘UB))*. (42) 

If we multiply both sides of Eq. (1421) by R n /{R u )r, with 
n = 1, 2, ..., d, and average over R we arrive at: 


t{n) 



(R n+d ~ k ) R {R k ) R 

(R d ) R (R n ) R 


(43) 


where 


t{n) 


{R n T(R)) R 

(R n )R 


(44) 


From Eqs. CM) and (1441) we see that the mean cluster 
size can be obtained from S = t( 0). 

To solve Eq. (l43l) , we note that for large d the binomial 
coefficient is strongly peaked at k = d/2 and takes the 
asymptotic form 




2 d +i 


2 k 1 \ 

( 45 ) 
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where g(x,cr) = exp(— x 2 / 2a 2 )/\/2t: a 2 is the Gaussian 
function. Provided that the radii distribution is bounded, 
the binomial coefhcent dominates the fc-dependence of 
the kernel. To see this, let us consider the m-th moment 
(■ R m )n = Rm fo d yp(y)y m i where y = R/R M ■ For large 
m the main contribution to the integral comes from y 
close to 1. Thus we make the quite general assumption 
that for y —> 1 , the radii distribution behaves as p(y) oc 
(1 — y ) a ~ 1 1 with a > 0. Setting t = m( 1 — y) for large 
m, we find 



d d 


ora /*ra 

{R m ) R oc —^ / dtt a ~ l {\ - t/m) r ' 

m J o 


o 

11 


M 

m a 


dtt { 


a ~ 1 e~ t = 


R 


M 

m a 


T(a), (46) 


so that for large k the term (R n+d ~ k ) R (R k ) R in Eq. fl3l) 
is proportional to R///~ d /[(n+d— k)k] a , which has a much 
weaker fc-dependence than Eq. (l45l) . Next, we introduce 
s = 2 n/d and s' = 2 k/d, which we treat as continuous 
variables for d —> oo, and we replace the sum over k 
by an integral over s': J 2 k =o “^ f Jo ds '• ^ we denote 
t(s) = t(ds/ 2) and t(s') = t(ds'/ 2), Eq. (TT7T1) becomes: 


t(s) = 1 + 2 d rj / ds 
Jo 


9 V lj Vd 

X (R d[ 1 +{s - s ' )/ 2 ] )R(R ds,/ 2 )R ^ r/) 


{R d ) R {R ds / 2 ) r 


(47) 


Since g(s' — 1, 1/y/d) —> S(s' — 1 ) for d —> oo, the above 
expression reduces to: 


t(s) = 1 + r ]2 


d (R d{ 1 +s ) / 2 )R( R d/2 )R 


(R d ) R (R d °/ 2 ) r 
from which we obtain the mean cluster size: 

S = t( 0) = 1 + r\2 d ^fp?” t(l). 


K 1), (48) 


(R d ) 


R 


(49) 


Setting s = 1 in Eq. C51) . we find t(l) = (1 — 2 d rj) 1 , so 
that we arrive finally at: 


i (R d/2 ) R 

1 — 2 d r] ( R d ) R ’ 


(50) 


which, as found for the case of discrete distributions, di¬ 
verges at 


b-^c = ^, (51) 

independently of the particular form of the bounded dis¬ 
tribution p{R). 

Using Eq. (1451) and considering the weak dependence 
of the moments of R, we readily obtain the large dimen- 


FIG. 4. (Color online) (a) Dimensional dependence of the 
percolation threshold g c in units of the asymptotic value 2~ d 
obtained from a numerical solution of Eq. (ITHl for rectangu¬ 
lar (upper curve), semicircular (middle curve), and triangular 
(lower curve) distributions of the sphere radii, (b) Corre¬ 
sponding critical average connectivity per particle z c (solid 
curves). Dashed lines are the asymptotic results for d 1 : 
z c = 4/d (rectangular distribution), z c = 32/(y^d 3 ^ 2 ) (semi¬ 
circular distribution), and z c = 32/d 2 (triangular distribu¬ 
tion). 


sional limit of the critical average connectivity per parti¬ 
cle: 


Zc = Pc^d{(R + R') d )R,R' 


( R d/2 )l 

(R d )R 


V- fd\ (R k ) R (R d ~ k )R 

win 

(52) 


from which we see that z c < 1 for any bounded distri¬ 
bution of the radii. Note that from Eq. (l52l) we recover 
z c = xm when p(R) is given by Eq. (I5U1) . 

We complete this section by showing how the per¬ 
colation threshold obtained from Eqs. flU and m 
evolves towards the asymptotic value r) c = 2 d as d 
increases. Toward that end, we consider radii dis¬ 
tributions of rectangular, semicircular, and triangular 
shapes, given respectively by p{R) = 1 /Rm, p{R) = 
A^(Rm/2 ) 2 + (R - Rm/^/tt, and p(R) = 2{R M - 
R)/R 2 m , for R < Rm and zero otherwise. We calculate 
T] c from the smallest pole of S = f (0) obtained from a nu¬ 
merical solution of Eq. (H51) . The resulting thresholds are 
very close to 2 ~ d for all d considered, and they approach 
the asymptotic limit from below, as shown in Fig. Oja). 
For the same cases of Fig. |4ja), we have calculated also 
the d-dependence of z c , shown in Fig. 0Kb) by solid lines, 
which we compare with the asymptotic limits (dashed 
lines) z c = 4 /d, z c = 32/( v / 7t d 3 / 2 ), and z c = 32/d 2 ob¬ 
tained from Eq. (1521) for rectangular, semicircular, and 
triangular distributions of the radii, respectively. 















VI. THE CASE OF UNBOUNDED 
DISTRIBUTION OF THE RADII 


Having established that r/ c is universal as d —> oo for 
bounded (and independent of d ) distributions of the radii, 
it is natural to ask if universality holds true also when 
p(R) is unbounded. Although we have shown the ir¬ 
relevance of closed loops limited to the case of bonded 
distributions, we shall nevertheless assume that n-cycle 
coefficients are negligible also for unbounded p(R), and 
that graph components have a tree-like structure. Let 
us consider the specific case of a lognormal distribution 
function: 


p(R) 


1 

yfTnoR 


exp 


ln(R/R 0 ) 
2 a 2 


(53) 


where R £ [0, oo), Rq is the median radius, and cr is the 
standard deviation of ln(i?). Equation (1531) is an interest¬ 
ing case-study, as the resulting ry c and z c for asymptot¬ 
ically large d can be found analytically. Using the k -th 
moment ( R k )n = Rq exp(cr 2 fc 2 /2), Eq. (H51) becomes: 


t(n) = 1 + V (fc) e^ 2 ^ n+d - k) 2 +k 2 - n 2 - d 2 ] t{k), (54) 
from which we express the mean cluster size as: 

S' = f (0) = 1 + 77 V ( f ) e a 2 k(k - d) t[k). (55) 

k—0 ' ' 

For sufficiently large d, the only nonvanishing terms of 
the summation are those with k = 0 and k = d, so that: 

S = 1 + r][S + t(d)}, (56) 

where from Eq. (1541) t(d) is given by: 

t(d) = l + V p^e° 2 ( d -V 2 t(k). (57) 


For d —> oo, t(d) tends asymptotically to t(d) = 1 + 
T/e * 7 d t( 0 ), as the term with k = 0 dominates the sum 
over k in Eq. m- We thus find that the mean cluster 
size, Eq. (|56]l . reduces to: 


S = 


1 + 77 

1 — 77 — rj 2 e a2d2 ’ 


(58) 


which diverges at the asymptotical critical value, 


_ 1 -2 .2 
r/^r) c = e ^ a 


(59) 




FIG. 5. (Color online) (a) Dimensional dependence of the 
percolation threshold r/ c in units of the asymptotic value 
exp(— a 2 d 2 /2) for a lognormal distribution of the radii ob¬ 
tained from a numerical solution of Eq. (f54|) : cr = 0.25, 0.3, 
0.4, and 0.5 from the lowermost to the uppermost curves, (b) 
Critical average connectivity per particle z c for a = 0.25, 0.3, 
0.4, and 0.5 from the uppermost to the lowermost curves. All 
curves tend to Zc/rjc —> 2 as d —> 00 . 


where we have again used the fact that for large d only 
the terms k = 0 and k = d contribute to the summation. 

As evidenced in Eq. (15^1) . the percolation threshold 
for infinite dimensions is no longer universal, as it de¬ 
pends on the parameter cr of the log-normal distribu¬ 
tion. Interestingly, from Eq. (1591) we also see that r] c can 
be smaller than the critical threshold of monodisperse 
spheres ( rj c = 2~ d ), contrary to what is expected in finite 
dimensions [24[ . We note that a critical threshold smaller 
than the monodisperse sphere limit in large dimensions 
has been found also for the case of radii distributions 
with d-dependent weights [ 25 ], |26| . 

To verify the accuracy of Eq. (1591) , we compare it with 
the threshold obtained by solving numerically Eq. m- 
As d increases, the asymptotic limit rj c = e~^ a d is 
reached more rapidly when cr is larger, as shown in 
Fig. EKa). From inspection of Eq. (l55l) we see that this 
behavior is due to the competition between e a fc ( fc-d ) and 
the maximum value ~ 2 d of the binomial coefficient at 
k = d/ 2 : the latter is suppressed by the exponential func¬ 
tion when d > 41n(2)/cr 2 . From numerical calculation of 
z c , shown in Fig. [51(b) for the same a values of Fig.[5](a), 
we see that also the asymptotic formula for z c , Eq. m, 
is verified. 


VII. LOWER BOUND ON THE PERCOLATION 
THRESHOLD 


The corresponding critical coordination number is 


z c 


V- fd\ ( R k ) R (R d ~ k ) R 
<*>» 



2??c, 


(60) 


Having established that Eqs. m and (1991) give asymp¬ 
totic limits of the critical threshold rj c as d —> 00 , we 
show now that the same equations provide also a lower 
bound on r] c for any dimensionality. Toward that end, we 
take the pair-connectedness function Pij (r) considered in 
Sec. IIV Cl and we extend to the polydisperse sphere case 
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the inequality formulated in Ref. f43| : 

p ij ( Y ) < fair) +pJ2 Xl I dr 'M\ r ~ r 'l) P b( r '); ( 61 ) 

i d 

where fij(r ) is the connectedness function given in 
Eq. ©. The above expression applies to any dimension¬ 
ality, and following Ref. jUJ, where Eq. m has been 
used for the monodisperse sphere case, it enable us to 
find a lower bound on the percolation threshold. To see 
this, we take the volume integral of Eq. (ED, 

Pij^V^+p^XiV^Pij, (62) 

i 

where = J drfij(r ) = Qd(Ri + Rj) d , and we use 
Eqs. (E5l) to find: 

R < 1 +pY,x j V£T j = 1 + E( 63 ) 

3 3 

which together with Eq. (f26ll gives an upper bound for 
the mean cluster size: 

s = Y, x iTi<Y, x * T » (64) 

i i 

where Ti is the solution of Eq. m- From the inequal¬ 
ity of Eq. (l6dl) . we see that the value of r] such that 
]Tb XiTi diverges identifies a lower bound on the perco¬ 
lation threshold for any d. The solid lines plotted in 
Figs. [3])a)-[5Ka) represent thus lower bounds on r) c for the 
different radii distribution functions considered in this 
work. As d increases, these lower bounds tend asymptot¬ 
ically to the infinite dimensional limit 2~ d for bounded 
radii distributions and to exp(— cr 2 d 2 /2) for lognormal 
radii distributions. Finally, we note that Eq. (1641) im¬ 
plies also that the values of z c shown in Figs. 0b)-[5](b) 
are lower bounds on the critical connectivity for any di¬ 
mensionality. 

VIII. SUMMARY AND DISCUSSION 

We have considered random dispersions of penetrable 
d-dimensional spheres with distributed radii in terms of 
weighted random geometric graphs, where nodes repre¬ 
sent sphere centers and edges connect nodes of overlap¬ 
ping spheres with probability weighted by the sphere 
radii. For bounded distribution of the radii, we have 
shown that closed loops of connected spheres can be ne¬ 
glected in the limit d —> oo and that graph components 
have thus tree-like structure. Analysis of the mean clus¬ 
ter size reveals that the asymptotic percolation threshold 
is universal and coincides with the threshold rj c = 2~ d 
found for the case of monodisperse spheres in high di¬ 
mensions. This result confirms and extends a previous 
finding on the percolation of d —> oo spheres with two 
different radii [25l. [2(1 . Furthermore, we show that the 


asymptotic critical connectivity per particle z c , though 
dependent on the shape of the radii distribution function, 
is less than unity and approaches z c —> 1 for spheres of 
identical radii. 

We have also studied critical connectivity for spheres 
with radii distributed according to a ^-independent log¬ 
normal function, which is a treatable example of un¬ 
bounded distribution. Assuming that clusters have a 
tree-like structure, we find that the percolation threshold 
t/ c depends on the shape of the log-normal distribution 
and, interestingly, that r] c for d —> oo can be smaller that 
the threshold for monodisperse spheres, in contrast to 
what is expected at finite dimensions |24j . 

Before concluding, let us speculate on the percolation 
threshold in homogeneous fluids of polydisperse spheres 
with impenetrable cores (cherry-pit model |3(). In finite 
dimensions, correlations between the cores preclude writ¬ 
ing the multi-degree distribution as a product of Pois¬ 
son distributions, as done in Sec. IIV A1 because the N- 
particle distribution function gjy( ri, r 2 ,..., r^r) depends 
on the relative positions of the core centers |32|. However, 
in the limit of infinite dimensions and for small densities, 
< 7 ,/v(ri, r 2 ,..., r^r) asymptotically factorizes into a prod¬ 
uct of 0-functions that are unity for pair distances beyond 
the hard-core diameter |44j]. The multi-degree distribu¬ 
tion for d —> oo can thus still be written as a product of 
Poisson distributions, with the average number of con¬ 
tacts unaltered by the presence of the hard-cores if the 
penetrable shells are non-vanishing. With the same rea¬ 
soning, closed loops are expected to be negligible and 
graphs are still dominated by tree-like components. For 
non-zero penetrable shells, therefore we expect the same 
asymptotic results for r] c as those obtained for the case 
of penetrable hyperspheres. 
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Appendix A: Irrelevance of (cd)^ for d — > oo 

In this appendix, we show that when the radii distri¬ 
bution is independent of d and bounded [that is, when 
p(R) = 0 for any R > Rm, with Rm < oo], the n-cycle 
coefficient for polydisperse spheres, defined in Eqs. 0- 
©. vanishes for d — > oo. 

Since Rm is the maximum radius of the distribution, 
the connectedness functions in the integrand of Eq. 0 
are such that fij(r ) < f(r) = 9{2Rm — r) for any i and 
j. We can thus write: 

C nl.,i n < [ dr {n) f( |ri - r 2 |)/(|r 2 - r 3 |) • • • /(|r„ - ri|), 

(Al) 
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which, when substituted in Eq. 0, gives: 


noting a general ra-th moment as we obtain: 


(c d ) {n) < C, 


(n) 


vv 

v v p 


n -1 
ex 




( n) 


X- 


(A2) 


where V ex = 2 d CldR%, and c£ l ' > is the n-cycle coefficient 
for identical radii given in Eq. ©• Next, we perforin the 
integrations over r 1; ..., r n in Eq. © to find: 




= vnr 


d 

kn— 1 


'£(0- £ 

k i=0 V 7 fc„_i=0 

x {R kl ) R (R d - kl+k2 ) R (R d ~ k2+k3 ) R ■ ■ ■ 

■ ■ ■ (A4) 


Following Sec. IV B1 we approximate for large d the bino¬ 
mial coefficients by Gaussian functions centered at d/2, 
and we replace the sums by integrals, so that for d —> oo, 
E* J= o (j$) 2d /o dajS{aj - 1), where a, = 2 kj/d and 
j = 1 ,,n— 1. Equation (1X1 reduces in this way to: 


<v£„ -»• v^-^-^r^KrX- 2 3 4 5 6 7 8 9 10 11 ^ 

(A5) 

so that Eq. (IA2I) becomes: 


V, 


(n) 


n —1 n —1 

=^ n ^ ,ij+i =^r 1 n x + R ^° 

i=i i=i 


= ^r i n 

1=1 


E (/)<'<« 

fc,=o 


(A3) 


where in the last equality we have expanded the binomial 
powers. In performing the average over R .^,..., Ri n , we 
must group the contributions with equal radius variables 
and average them independently of the other radii. De- 


(q)<"> < 4"M n) , 


(A6) 


where 


(n) 

Xd = 


i? 


(n—l)d 

■M 


(R d /*)URT R 2 ' 


(A7) 


For continuous radii distributions, we assume that 
p(R) oc (Am — A)“ -1 with a > 0 for R —1 Rm, as done 
in Sec. IV 61 Using Eq. (H>)) we thus find oc d na . For 
discrete distributions as in Eq. (Id()l) , it is straightforward 
to show from Eq. (IA7D that Xd' > 1 /(*m)” for large 

(n) 

d. We have thus arrived at the result that \d increases 
with d at most as a power-law, leading to ( c d 1 0 as 
d —1 oo, due to the exponential vanishing of . 
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